Light Cone Sum Rules for the 7r°7*7 Form Factor Revisited 
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We provide a theoretical update of the calculations of the 7r°7*7 form factor in the LCSR frame- 
work, including up to six polynomials in the conformal expansion of the pion distribution amplitude 
and taking into account twist-six corrections related to the photon emission at large distances. The 
results are compared with the calculations of the B —> -kIv decay and pion electromagnetic form 
factors in the same framework. Our conclusion is that the recent BaBar measurements of the 7r°7*7 
form factor at large momentum transfers [lj are consistent with QCD, although they do suggest 
that the pion DA may have more structure than usually assumed. 
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I. INTRODUCTION 

Despite a solid theory background [2-4], phenomeno- 
logical success of QCD in exclusive reactions has been 
rather modest. A problem is that, since the quarks carry 
only some fractions of hadron momenta, virtualities of 
the internal lines of the hard subprocess appear to be es- 
sentially smaller than Q 2 , the nominal momentum trans- 
fer to the hadron. As a result, at accessible Q 2 , the 
bulk part of the hard QCD contribution comes from the 
regions where the "hard" virtualities are much smaller 
than the typical hadronic scale of 1 GeV 2 043] • Accord- 
ing to the factorization principle, contributions from such 
regions have to be subtracted from the hard coefficient 
function and included separately as additive "soft" or 
"end-point" nonperturbative contributions. The stan- 
dard power counting suggests that "soft" terms are sup- 
pressed by extra powers of l/Q 2 . However, they do not 
involve small coefficients ~ a s (Q)/ir which are endemic 
for factorizable QCD contributions based on hard gluon 
exchanges. As the result, the onset of the perturbativc 
regime may be postponed to very large momentum trans- 
fers. 
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least one virtual photon plays a very special role as it 
is the simplest hard exclusive process where the above 
mentioned difficulties are absent or, at least, moderated. 
There is only one hadron (pion) involved, and the large 
Q 2 behavior of this form factor is determined [4j by the 
operator product expansion (OPE) of the product of two 
electromagnetic currents near the light cone which is very 
well studied in the context of inclusive deep-inelastic scat- 
tering (DIS). In this case the leading contribution to the 
hard coefficient function is of order one (i.e. not sup- 
pressed) as no gluon exchanges are involved, and at the 
same time "soft" (end-point) contributions either do not 
exist — for the case of two virtual photons — or are likely 
to be suppressed, if one photon is real. These features 
make the pion transition form factor an ideal place to 
test the QCD factorization approach and determine the 
pion distribution amplitude (DA) which can then be used 



to describe other exclusive hard reactions. This task is 
as important as ever, the most high-profile application 
being at present to exclusive weak B-decays, B — > nivi, 
B — > 7T7T etc. which are the main source of precision in- 
formation on quark flavor mixing parameters in the Stan- 
dard Model. These are the aim of an extensive experi- 
mental study: It addresses the question whether there 
is New Physics in flavor-changing processes and where it 
manifests itself. 

Whereas the case of two virtual photons offers crucial 
simplifications for the theory, the transition form fac- 
tors with one real and one virtual photon are much eas- 
ier to study experimentally. They can be measured for 
space-like momentum transfers in the process e + e~ — > 
e + e~7r°, 77, . . . and in e + e~ — > 777 ,rj,... for time-like 
ones. Till 1995 only the CELLO data [|[ were available 
at relatively low, space-like momentum transfer: Q 2 < 
2.5 GeV 2 for 7r° (77)7*7 an d somewhat higher for V7*7- 
The covered region was extended to Q 2 ~ 9— 15 GeV 2 by 
CLEO collaboration 0] which allowed, for the first time, 
a quantitative comparison with the perturbative QCD. 
The CLEO data appeared to be consistent with the pre- 
dicted scaling behavior ~ l/Q 2 setting in for momentum 
transfers of the order of a few GeV 2 and also suggested 
that the pion DA is somewhat broader compared to its 
asymptotic shape at large scales, which was, again, ex- 
pected based on the corresponding studies using QCD 
sum rules. More recently, BaBar reported [l(| a mea- 
surement of the time-like transition form factors 777*7 
and ^'7*7 at very large q 2 — 112 GeV 2 . No significant 
tension with the theory was observed, although predic- 
tions in the time-like region are generally more difficult. 

The situation changed in 2009 when the BaBar collab- 
oration presented [l| the measurement of the 7r°7*7 form 
factor up to photon virtualities of the order of 40 GeV 2 . 
These new data created considerable excitement in the 
theory community as they do not show the expected scal- 
ing behavior. The most popular explanation so far has 
been [ll|, that the pion DA has an unexpected "flat" 
shape and does not vanish at the end points. This, on 
one hand, triggered speculations on the breakdown of 
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QCD factorization [l3| and, on the other hand, claims 
that the BaBar data are in contradiction with light-cone 
sum rules (LCSRs) and with the common wisdom on the 
lowest moments of the pion DA [H-[l6|. The aim of this 
work is to reexamine these claims by making an updated 
analysis of the tt°j*j form factor within the LCSR ap- 
proach including up to six polynomials in the conformal 
expansion of the pion distribution amplitude and taking 
into account twist-six corrections. The photon emission 
at large distances is discussed in detail. The results are 
compared with the calculations of pion electromagnetic 
form factor and B — > iriv decay in the same framework. 
Our conclusion is that the recent BaBar measurements 
PJ are consistent with QCD and with the bulk of the 
available information on the pion distribution amplitude. 
In particular we argue that the "flat" DA [TO, Qj] is not 
warranted and in fact no conclusion on the end-point 
behavior of the DA can be inferred on the basis of the 
existing experimental data. 

The presentation is organized as follows. Sect. 2 con- 
tains a concise review of the QCD (collinear) factoriza- 
tion approach to the 7r°7*7 form factor and the existing 
information on the pion DA. The LCSR approach is mo- 
tivated and explained in detail in Sect. 3. In this work 
we go beyond the existing analysis 0, [l7- 22] in two as- 
pects. First, we calculate a new, twist-six contribution to 
LCSRs which proves to be sizeable. This contribution is 
related to photon emission from large distances for which 
we also derive the leading-order perturbative expression. 
Second, we extend the existing formalism to allow for the 
contributions of higher-order Gegenbauer polynomials, 
which allows one to consider DAs of arbitrary shape and 
also address the question of convergence of the Gegen- 
bauer expansion which generated some confusion. The 
second task was already addressed in Ref. [14j |. but our 
expressions do not agree, unfortunately. Sect. 4 contains 
the numerical study of the LCSRs to the NLO accuracy. 
We consider various uncertainties of the method in some 
detail and provide error estimates for our predictions. 
The final Sect. 5 is reserved for a summary and conclu- 
sions. 



II. QCD FACTORIZATION AND PION DA 



The form factor F 1 * 1 *^o{q\ , q%) describing the pion 
transition in two (in general virtual) photons can be de- 
fined by the matrix element of the product of two elec- 
tromagnetic currents 



d 4 ye^y(n\p)\T{f™(y)jr(0)}\0) 



= «e 2 e Aira/3 (?"^F 7 . 7 ,^ 7r o(gi,(j|) (1) 

where 

JT(y) = e«u(y)7/Xy) + e d d(yhf,d(y) + 

p is the pion momentum and q2 = qi + P- We will con- 
sider the space-like form factor, in which case the both 



photon virtualities are negative. The experimentally rel- 
evant situation is when one virtuality is large and the 
second one small (or zero). For definiteness we take 
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assuming that q 2 -C Q 2 ■ Most of the equations are 
written for q 2 = and we use a shorthand notation 

F^^oiQ 2 ) = F^^ n a(q 2 = -Q 2 ,q 2 = 0). 

In general, a power-like falloff of the form factor 
F 1 * 1 ^ yK o{Q 2 ) in the large-Q 2 limit can be generated by 
the three different possibilities of the large- momentum 
flow as indicated schematically in Fig. [T] 23] . The first 
possibility, Fig. QJi, corresponds to the hard subgraph 
that includes both photon vertices. This is the domi- 
nant contribution that starts at order ~ 1/Q 2 - For zero 
(or small, q 2 < Aq CD ) virtuality of the second pho- 
ton there exists another possibility shown in Fig. [Tp: 
In this case the low-virtual photon is emitted at large 
distances and the large momentum flows through a sub- 
graph corresponding to hard gluon exchange between the 
quarks. The power counting for this contribution shows 
that it is at most 0(1/Q 4 ), i.e. subleading compared 
to the first regime. Finally, the third possible regime 
shown in Fig.[T]; corresponds to the Feynman mechanism, 
i.e. the possibility that the quark interacting with the 
hard photon carries almost all the momentum whereas 
the quark spectator is soft. This contribution can intu- 
itively be thought of as an overlap of nonperturbativc 
wave functions describing the initial (photon) and final 
(pion) states. In perturbation theory, this contribution 
also scales as 0(1/Q 4 ) in the large-Q 2 limit. 

The contribution in Fig. QJ, by construction involves 
a time-ordered product of two electromagnetic currents 
at small light-cone separations. Hence it can be studied 
using Wilson operator product expansion. The leading 
contribution 0(1/Q 2 ) to the form factor corresponds to 
the contribution of the leading twist-two operators and 
can be written in the factorized form 



^7*7->7r°(Q 2 ) 



dxT H (x,Q 2 , n,a s (n))4>n(x, m) j 

(3) 

where ff (x,//), the pion distribution amplitude at the 
scale [J., is defined by the matrix element of the nonlocal 
quark-antiquark operator stretched along the light-like 
direction n^, n 2 = 0: 

(0|g(0)[0, cm] rt l5 q(an)\n+ (p)> = 



= iUp-n / dxe- lxap - n ^(x,fi). (4) 
Jo 

Here and below qifi^q = (l/\/2) [1/^7511 — d'dr/^d]. To 
this accuracy (leading twist) all gluon attachments to 
the hard subgraph (cf. Fig. [TJl) can be absorbed in the 
path-ordered gauge link (Wilson line): 

[0,cm] = Pexp|— ig J du n^A^ (un) | . (5) 
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The normalization is such that 



1 l9(0)7.75g(0)|7r°(p)) = if vPv , 



dx <^>7r(a;, n) 



(6) 

where f„ ~ 131 MeV is the usual pion decay constant. 

The coefficient function in Q is known in the MS 
scheme to the next-to- leading order (NLO) in the strong 
coupling [24T4261 ] . Taking into account the symmetry of 
the pion DA (^^(x) — </>?,- (1 — x) it can be written as 



T, 



NLO 



H 



xQ 



Q 2 i 



l + C F - 
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lna; In 



Q 2 



:ln: 



2(1 -x) 



(7) 



Symmetry properties of the renormalization-group 
(RG) equation which governs the scale dependence of the 
pion DA 0, 3] suggest the series expansion of the DA in 
Gegenbauer polynomials 



71=0 



ip n (x) = Qx(l- x)C 3 J 2 (2x-l) 



(8) 



The first coefficient ao(/i) = 1 is fixed by the normal- 
ization condition whereas the remaining ones, a„(/io) for 
n = 2, 4, . . ., have to be determined by some nonpertur- 
bative method (or taken from experiment). 

To leading order (LO) the Gegenbauer coefficients are 
renormalized multiplicatively whereas to the NLO accu- 
racy the mi xing pattern becomes more complicated. One 
obtains l27H32f 



n-2 



4tt 



^a k ((io)El (fx,» )d*(ix,ii ), (9) 



fe=0 



Explicit expressions for the RG factors E^ LO (fj,, ^o) and 
the off-diagonal mixing coefficients d^(/x,^ ) i n the MS 
scheme are collected in Appendix A. 

The NNLO calculations of the transition pion form fac- 
tor exist in the so-called conformal scheme CS 33, 34] but 
they cannot be converted to MS lacking the full NNLO 
result for the trace anomaly term, which is so far not 



available. An extensive discussion of scheme dependence 
can be found in Refs. [31TI33} . 

The existing information on the pion DA comes from 
QCD sum rules, lattice calculations and light-cone sum 
rules. The first nontrivial Gegenbauer coefficient a 2 is 
related to the second moment of the DA 



(£ 2 ) = / dx(2x- lf^(x). 
Jo 



a 2 



(10) 

and can be evaluated as a matrix element of the local op- 
erator with two derivatives between vacuum and the pion 
state. There exists overwhelming evidence that this co- 
efficient is positive, meaning that the pion DA is broader 
than its asymptotic expression </>^ s = 6a; (1 — x), see Ta- 
ble u 

Such calculations where pioneered in 1981 by Chernyak 
and Zhitnitsky 35] who derived the corresponding sum 
rule and obtained a 2 ~ 0.5 at the scale of order ^l 2 = 
1 — 1.5 GeV 2 . Extrapolating this number to a very 
low scale fi 2 = 0.25 GeV 2 and adding a simplifying 
assumption that higher-order coefficients 04,... vanish, 
they have formulated a simple model for the low-energy 
pion DA, which has become known as the Chernyak- 
Zhitnitsky (CZ) model: 



cj)^ z (x) = 30x{l - x)(2x-iy 



6x(l — x) 



l+ l -CT(2x 



(11) 



This model corresponds to a very asymmetric momentum 
fraction distribution which vanishes at the point where 
the pion momentum is shared equally between the quark 
and the antiquark <^ z (x = 1/2) =0. The striking differ- 
ence of the CZ model and the asymptotic DA has been 
fuelling an extensive and sometimes heated discussion for 
many years. 

Newer estimates of a 2 following the CZ approach yield 
a somewhat smaller value [33, [3^] 02(1 GeV) ~ 0.3, the 
difference being due to a combination of several factors: 
writing the sum rules for a 2 directly instead of the second 
moment (£ 2 ), taking into account the NLO corrections 
and using slightly different values of the parameters. 

Light-cone sum rules J49-[51J are a modification of the 
general SVZ approach (52|, in which the pion DA serves 
as the main input in calculations of form factors (or 
hadron matrix elements) . The Gegenbauer coefficient a 2 
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Method 


fj, = 1 GeV 


H = 2 GeV 


Reference 


LO QCDSR, CZ model 


0.56 


0.39 


[35, 36] 


QCDSR 


n 26+ 21 

u - zo -0.09 




[37] 


QCDSR 


0.28 ±0.08 


0.19 ±0.06 


m 


QCDSR, NLC 


0.19 ±0.06 


0.13 ±0.04 


[19, 39, 40] 


F n ryy* j LCSR 


0.19 ±0.05 


0.12 ±0.03 Gu = 2.4) 


[18] 


F n ryy* j LCSR, 


0.32 


0.21 = 2.4) 


[20] 


■^Tvyy* i LCSR., R 


0.44 


0.31 


[41] 


J^7T77* 5 LCSR, R 


0.27 


0.19 


[22] 


FT\LCSR 


0.24 ±0.14 ±0.08 


0.17 ±0.10 ±0.05 


[42, 43] 


-F^ m ,LCSR, R 


0.20 ±0.03 


0.14 ±0.02 


[44] 


Fb-^ttIv, LCSR 


0.19 ±0.19 


0.13 ±0.13 


[45] 


Fb^ttIv, LCSR 


0.16 


0.11 


[46] 


LQCD, N f = 2, CW 


0.289 ±0.166 


0.201 ±0.114 


QCDSF/UKQCD [47] 


LQCD, N f = 2+1, DWF 


0.334 ±0.129 


0.233 ±0.088 


RBS/UKQCD [48] 



TABLE I: The Gegenbauer moment a2(/i 2 ). The CZ model involves ai = 2/3 at the low scale fi = 500 MeV; for the discussion of 
the extrapolation to higher scales, see Ref. [201] . The abbreviations stand for: QCDSR: QCD sum rules; NLC: non-local condensates; 
LCSR: light-cone sum rules; R: renormalon model for twist-4 corrections; LQCD: lattice calculation; Nf = 2(±1): calculation using 
Nf = 2(+l) dynamical quarks; CW: non-perturbatively O(a) improved Clover-Wilson fermion action; DWF: domain wall fermions. 
For convenience we present the results for two scales, fj, — 1 GeV and /j, = 2 GeV, the relation is calculated in NLO. 



is not calculated directly, but is extracted from the com- 
parison of the LCSR calculations with the experimental 
data. Note that in the case of the pion transition form 
factor these fits are based on CLEO data [9] only. The 
results for a-i are consistent with the direct calculations, 
see Table [D 

Finally, two independent lattice calculations of ai are 
now available [13, EH • The largest part of the uncertainty 
in these results is due to the chiral extrapolation. Overall, 
the results in Table U show a consistent picture 

4 LO C« 2 = 1 GeV 2 ) = 0.25 ± 0.10 (12) 

and one can expect that the accuracy will increase in near 
future when lattice calculations with physical pion mass 
become available. 

Very little, unfortunately, is known about the next 
Gegenbauer coefficient, 04. The LCSR fits of heavy me- 
son decay form factors indicate a small positive value, 
04 ~ 0.04 J46|, whereas the similar approach applied 
to pion transition form factor (CLEO data only) favors 
small negative values [2(| [22], |4l[ . Lattice calculations of 
this coefficient suffer from large (lattice) artifacts in the 
operator renormalization and are not feasible at present. 

The calculation of within the QCD sum rule ap- 
proach has been attempted in the so-called nonlocal con- 
densate model (NLC) [HIIIlII which involves a rcsum- 
mation of a tower of condensates of a certain type. This 
approach leads to a sizeable negative value 04 ~ —0.1 
which is included in the so-called Bakulev-Mikhailov- 
Stefanis (BMS) model of the pion DA. A large negative 
value for 04 in the NLC approach can be traced back 
to the basic feature of this model that nonperturbative 



corrections to the DA get smeared over a finite interval 
of momentum fractions Ax ~ A 2 /(2M 2 ) ~ 0.2 where 
A 2 = (qD 2 q) / (qq) ~ 0.4 GeV 2 is the average virtuality of 
quarks in QCD vacuum and M 2 ~ 1 GeV 2 is the Borcl 
parameter. To our opinion, appearance of Air is an arti- 
fact of the NLC model: contributions of this type would 
produce "bumps" at large values of Bjorken variable in 
quark parton distributions in the nucleon [53] (which are 
absent) and also a finite smearing proves not sufficient 
to cure the QCD sum rules for heavy-to-light decay form 
factors (which is the reason why this technique was 
eventually abandoned and replaced by LCSRs). It seems 
much more natural to assume that the nonperturbative 
contributions get smeared over the whole interval of mo- 
mentum fractions < x < 1. A possible mechanism for 
such "complete" smearing is considered on the example 
of a photon DA in Appendix B in Ref. [Hj]. We believe, 
therefore, that the NLC-model-based predictions for 04 
have to be viewed with caution. 

Last but not least, we mention the LCSR calculation 
[Hol ] for the pion DA in the middle point: 

<p w (x = 0.5, ^ 2 = 1 GeV 2 ) = 1.2 ±0.3, (13) 

which is only available result beyond the Gegenbauer ex- 
pansion. This result excludes a large "dip" in the pion 
DA in the center region and thus contradicts the CZ and 
BMS models. It is, however, consistent with most of the 
parameterizations of the pion DA that are used in vast 
literature on B-decays. Smaller values of ^(x = 0.5) are 
also strongly disfavored by numerous LCSR calculations 
ofpion-hadron couplings j^jvjv, g-^DD' , 9 P ^, see e.g. 
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FIG. 2: Pion transition form factor in the MR model (fTl)) 
calculated using flat pion DA (solid red curve) and CZ-type 
DA with <i2 = 0.5 and higher Gegenbauer moments set to 
zero (dashed blue curve). The nonperturbative parameter 
a = 0.53 GeV 2 in both cases. The experimental data are 
from [H (full circles) and [U (open triangles). 



the data, see the dashed (blue) curve in the same fig- 
ure. Moreover, it is seen that the two DA models can 
hardly be discriminated at all unless precise data with 
Q 2 > 20 GeV 2 are available! 



easy to see why this happens. The fiat DA 
k (x) = 1 can be expanded in Gegenbauer polynomials 
as follows: 



It is 

/.fiat/ 



E 



,flat 



<Pk{x) , 



fe=0,2, 



flat _ 2(2fc + 3) 
k 3(fc + l)(fc + 2) 



(15) 

This equation has to be understood in the sense of dis- 
tributions: The equality holds when both sides are in- 
tegrated with a test function that is finite (or does not 
increase too fast) at the end points. 

In perturbation theory (LO) the form factor is propor- 
tional to the sum of the Gegenbauer coefficients 



It has been suggested (TT|) EH that the BaBar data [l| 
indicate an unusual "flat" DA (f> n (x) ~ 1 which does 
not vanish at the end points so that the Gegenbauer ex- 
pansion is not convergent (more precisely: not uniformly 
convergent at the end points). We believe that this con- 
clusion is not warranted and in fact no conclusion on the 
end-point behavior of the pion DA can be inferred on the 
basis of the existing experimental data. 

To explain our statement, let us examine the argu- 
mentation in [Tl[ more closely. It is based on the ele- 
gant model for the transition form factor suggested by 
Musatov and Radyushkin (MR) in an earlier work (23j . 
which is derived from the exact two-body (e.g. quark- 
antiquark) contribution in the noncovariant light-cone 
formalism of Brodsky and Lepage Q under certain sim- 
plifying assumptions on the shape of the pion wave func- 
tion: 



pMR 

7*7- 



.AQ 2 ) 



V2U 
3Q 2 



dx 



<t>Tt{x) 



1 — exp 



X' 



IV 

2x<t J 



The first term in the square brackets corresponds to the 
usual leading-order contribution to Fig. [TJi, whereas the 
second term is entirely a soft contribution of the type 
Fig- [It- Note that this second term is exponentially small 
in Q 2 for each finite quark momentum fraction, so it is 
absent in any order of the operator product expansion. 
Using Eq. with a = 0.53 GeV 2 [H| and the "flat" 
pion DA (f)„(x) = 1 allows one to describe the apparent 
scaling violation in the BaBar data, as illustrated by the 
solid curve in Fig. [5J 

The caveat with this argument (and a very similar ar- 
gumentation in Ref. 12]) is that flat DA is not necessary; 
in fact a CZ-type DA with a 2 = 0.5 



<t>ir{x) — 6x(l — x) 



l+ l -Cl /2 {2x~l) 



1 dx 
o x 



1 + a 2 + a 4 + 



For the flat DA this series diverges, which motivates in- 
troduction of a certain regulator, e.g. the soft correction 
given by second, exponential, term in Eq. (|14p or effec- 
tive quark mass in Ref. [l2|. Our observation is, however, 
that if such a regulator is included, the Gegenbauer ex- 
pansion for the form factor is converging very rapidly and 
contributions of higher-order terms in this series are in 
fact negligible. For illustration, consider the MR model 
with a flat DA for, say Q 2 = 20 GeV 2 , and check how 
much is being contributed by each successive Gegenbauer 
polynomial. One obtains 



Q 2 ^"o(Q 2 = 20) = 

V2fn 



V2fn 



• 3.56513 



2.724 + 0.649 + 0.162 + 0.028 



(16) 



yields a very similar (or even better) description of 



where the first term on the r.h.s. is the contribution of 
the asymptotic DA, the second term is due to af 31 etc. 
One sees that all contributions beyond n — 4 are tiny. 

Our conclusion is that a good description of the BaBar 
data [l[ achieved in (Til [l2[ is not due to an unusual 
end-point behavior of the DA, but rather to a (model 
dependent) large nonperturbative soft correction to the 
form factor. Such a large correction effectively suppresses 
contributions of higher order terms in the Gegenbauer ex- 
pansion and makes the question of the end-point behav- 
ior of the pion DA irrelevant. The problem is, therefore, 
whether such a large nonperturbative correction can be 
expected, and whether it can be estimated in a model- 
independent way. This is the question that we address 
in the next Section. 
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III. LIGHT CONE SUM RULES FOR THE 
PION-PHOTON TRANSITION 

A. Dispersion approach 

The technique that we adopt in what follows has been 
suggested originally by Khodjamirian in [17]. It is to 
calculate the pion transition form factor for two large 
virtualities, Q 2 and q 2 , using the OPE, and make the 
analytic continuation to the real photon limit q 2 = us- 
ing dispersion relations. In this way explicit evaluation 
of contributions of the type in Fig. [TJd,c is avoided (since 
they do not contribute for sufficiently large q 2 ) and ef- 
fectively replaced by certain assumptions on the physical 
spectral density in the g 2 -channel. 

The starting observation [l?} is that F 7 » 7 «^„.o (Q 2 , q 2 ) 
satisfies an unsubtracted dispersion relation in the vari- 
able q 2 for fixed Q 2 . Separating the contribution of the 
lowest-lying vector mesons p, u> one can write 



F^,^^o(Q , q 2 ) 



7 7 



+ 



V2f p F Tp ^o(Q 2 ) 



s + q l 



where sq is a certain effective threshold. Here, the p 
and uj contributions are combined in one resonance term 
assuming m p ~ m u and the zero-width approximation 
is adopted; f p ~ 200 MeV is the usual vector meson 
decay constant. Note that since there are no massless 
states, the real photon limit is recovered by the simple 
substitution q 2 -> in (fT7|) . 

On the other hand, the same form factor can be calcu- 
lated using QCD perturbation theory and the OPE. The 
QCD result satisfies a similar dispersion relation 



2 „2\ 



1 



ds ■ 



lmF^ n0 (Q 2 ,-s) 



(18) 

The basic assumption of the method is that the physical 
spectral density above the threshold so coincides with the 
QCD spectral density as given by the OPE: 



lmF rr ^o(Q 2 , - a ) = lmFQ c ° (Q 2 , 



(19) 



for s > sq. This is the usual approximation of local 
quark-hadron duality. 

We expect that the QCD result reproduces the "true" 
form factor for large values of q 2 . Equating the two rep- 
resentations (fT7 |) . (fl"8|) at q 2 — > — oo and subtracting the 
contributions of s > sq from both sides one obtains 

V2/ p F 7W o(Q 2 ) = - / °dsImFQ c » (Q 2 ,~ S ) .(20) 

71" JO 

This relation explains why sq is usually referred to as the 
interval of duality (in the vector channel). The (pertur- 
bative) QCD spectral density ImF®, ^ t0 (Q 2 , — s) is a 



smooth function and does not vanish at small s — >• 0. 
It is very different from the physical spectral density 
ImF 7 » 7 *^ 7r o (Q 2 , — s) ~ S(s — m 2 ). However, the integral 
of the QCD spectral density over a certain region of en- 
ergies coincides with the integral of the physical spectral 
density over the same region; in this sense the QCD de- 
scription of correlation functions in terms of quark and 
gluons is dual to the description in terms of hadronic 
states. 

In practical applications of this method one uses an ad- 
ditional trick, borrowed from QCD sum rules [52j, which 
allows one to reduce the sensitivity on the duality as- 
sumption in Eq. (|19[) and also suppress contributions of 
higher orders in the OPE. The idea is essentially to make 
the matching between the "true" and calculated form 
factor at a finite value q 2 ~ 1 — 2 GeV 2 instead of the 
q 2 — >• oo limit. This is done going over to the Borel rep- 
resentation l/(s + <7 2 ) — > exp[— s/M 2 ] the net effect being 
the appearance of an additional weight factor under the 
integral 



V2f p F r * p ^o(Q 2 



dse -{s-ml)/M> 



xImF Q ^ (g 2 ,-,s).(21) 



Varying the Borel parameter within a certain window, 
usually 1-2 GeV 2 one can test the sensitivity of the results 
to the particular model of the spectral density. 

With this refinement, substituting Eq. (I21|) in (fl7|) and 
using Eq. (H"9l) one obtains for q 2 — > (l7j| 

F^ (Q 2 ) = i / S °4lmFQ C 7 ?_ (Q 2 ,- S )eK-)/M 2 



7T_ ^ 
1 

+ - 

7T 



ds , 



^ImF Q ™ o(QVs). (22) 



This expression contains two nonperturbative parameters 
- vector meson mass m 2 and effective threshold sq ~ 
1.5 GeV 2 — as compared to the "pure" QCD calculation, 
and the premium is that one does not need to evaluate the 
contributions of Fig. [T}d,c explicitly: They are taken into 
account effectively via the nonperturbative modification 
of the spectral density. 

As an illustration, consider the leading twist QCD ex- 
pression at the leading order 



F 



QCD , n 2 2n _ \/2/tt f 1 dx(j>„(x) 
7*7*->7r° W l" ) ~ o 



o xQ 2 + xq 2 



(23) 



In this case the momentum fraction integral can easily 
be converted to the form of a dispersion relation by the 
change of variables x —> s = Q 2 x/x. The resulting LO 
and leading twist LCSR expression is [l?} 



F LCSR (n 2\ 

7*7 — >7T f> 1^ ) 



V2U 



° dx <j>n(x) c (xm 2 r -xQ 2 )/(xM 2 ) 



o X "V 
1 dx cfrx (x) 



(24) 
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where 



x 



so 



(25) 



Note that the modification of the perturbative expression 
only concerns the region x < x$ ~ so /Q 2 so this is a soft 
contribution in our classification. Since 4>tt(x) ~ x for 
x — > 0, this contribution (first term in (|24jl ) corresponds 
to a power correction of the order of so/Q 4 for Q 2 — » oo, 
in agreement with usual reasoning based on the power 
counting. 



B. NLO perturbative corrections 

The NLO perturbative corrections to the 7*7*71-° form 
factor for arbitrary photon virtualities were calculated in 
Refs. 



CpOi s 
H — - — t(x, w) 



pQCD (n 2 2x _ V2U f 1 dx (j>„(x) 

,<Z j ~ 3 J xQ 2 +xq 2 



2tt 



(26) 

The function t(x, w) where w = Q 2 / (Q 2 + q 2 ) is given in 
Eq. (5.2) in [25j . 



Following Ref. [181 ] we write the required imaginary 
part of Fj, J?^a {Q 2 , q 2 ) as the sum of terms correspond- 
ing to the expansion of the pion DA cj) w (x, (J,) in Gegen- 
bauer polynomials ©: 



1 



ilmF^ (Q 2 ,-s) 
V2U 



(27) 



p1°>(QV) + ^p«(Q 2 , S ;m) 



where the LO partial spectral density is proportional to 
the pion DA 



pi 0) (Q 2 ,s) = 



Q 2 



Q 2 4 

The NLO spectral density 



dx ip n (x) — Im 



Q 



t(x, w) 



(28) 



xQ 2 +xq 2 



(29) 



can be written in the following form: 



Pn\Q 2 ,s) 



2(Q 2 4 

+ 7i 0) 



+ 2(V(2) - V(2 + n))] + y - In 2 (~) - 7 i 0) In (-^) } ^„(x) 



+ (x — > x) 



(30) 



where, as above, x = Q 2 /(Q 2 + s), tp(x) = dlxiT(x)/dx 
and 7^ is related to the leading-order anomalous dimen- 
sion 7i 0) (|A.4|) as 



(o) 



2C F 7, 



(31) 



Our result is similar but does not agree with the corre- 
sponding expression in Ref. The difference is that 
the first term in the second line in (|30l) is not sym- 
metrized in (x — > x) and hence the whole expression 
in braces is not symmetric under this substitution. We 
have checked that the spectral densities in (|30l) repro- 
duce the corresponding expressions for n = 0, 2, 4 in [18j : 

Pn\Q 2 ,s) = —A n . We also checked by numerical inte- 
gration for n < 12 that the dispersion relation 



ip n (x)t(x,w) 

dx —5 — 

x(y z -|-x<7 z 



S + q- 



•pW(g 2 , S ) (32) 



is indeed satisfied. 

As noticed in [lij], the integrals appearing in Eq. (|3"0)) 
can be expanded in terms of ip n (x) with rational coeffi- 



cients: 



l du Mu)-Vn{x) ^ 



n 

G*<p k (x), (33) 



fe=0.2, 



f/it 



¥>„(ll) - </?„(x) 



u — X 



-3x + ]T fl* Vfc (a;).(34) 



fc=0,l, 



The matrices and il„ can easily be calculated using 
orthogonality relations for the Gegenbauer polynomials, 
e.g. 



K = A/jT 1 f 1 dxC 3 k /2 (2x-l) 
Jo 



; v?„(u) - y> n (x) 
aw h 3x 



u — x 



(35) 



where A4 = (3/2) (A + l)(Jfe + 2)/(2fc + 3). Explicit ex- 
pressions for n, fc < 12 are collected in App. B. 
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C. Twist-four corrections 

Twist-four corrections to the form factor, by definition, 
correspond to the contributions of twist-four operators in 
the OPE of the time-ordered product of the two electro- 
magnetic currents in Eq. ([1]). Such contributions are of 
order 1/Q 4 , but, as we will see later, are not the only 
ones that have to be taken into account to this accu- 
racy: Twist counting does not coincide in the present 
case with the counting of powers of large momentum, so 
they should not be mixed. 

Intuitively, twist-four effects can be thought of as due 
to quark transverse momentum (or virtuality) in the 
handbag diagram shown in Fig. [3^, and quark- antiquark- 
gluon components in the pion wave function, Fig. [5Jd. 
These two contributions are related by exact QCD equa- 
tions of motion [6lJ; hence they must be taken into ac- 
count simultaneously. The corresponding calculation was 



done in Ref. 1_7J . The result (for two virtual photons) can 
be written as 



F » . 

^7 7 



V2U 



1 dx Mg) 

i, Q 2 x + q 2 x 



dx- 



o ~~{Q 2 x + q 2 x) 2 



(36) 



and therefore 



>ir°(Q 2 ,-s) 



V2U 



4>tt(x) 1 dFir(x) 



Q 2 Q 2 ds 



( 37 ) 

with the usual substitution x = Q 2 /(s + Q 2 ). The first 
term on the r.h.s. of ([55]). (|3"7|) is the leading order twist- 
two contribution and F w (x) is given in terms of the twist- 
4 pion DAs 



F 7r (x,/i) = -4>i^{x) + dai da 2 
4 Jo Jo 



1 / x— x+o 



a 3 



"3 /_ 



(38) 



where a = {ai, 0:2,0:3}. The two-particle twist-4 DA 
<t>i;-K{x) is defined by the light-cone expansion y 2 —> of 
the bifocal quark-antiquark operator [38l . l6ll [63 j 



|9(0)7m75«(2/)K(p)) = 

ri . 2 

= ifvpJ dxe~ ix ™ ^{x) + \-4>^{x) 
Jo 1 16 



.(39) 



and the three-particle DAs $4 ;7r (a), $4 ;7r (a) correspond 
to the matrix elements 

(0\q(Qhn"f5gG al 3(vn)q(un)\Tr(p)) = 

= Pn(p a np -ppn a ) — U®4w(u,v;pn) + 
pn 

(0\q(0)^igG a p(vn)q(un)\it(p)) = 

= P^Panp - pm a ) — f^A-Au, v : pn) + ... ,(40) 
pn 

with the shorthand notation 

The Wilson lines in the definitions of the nonlocal op- 
erators in (f39|) . (|40]) are not shown for brevity. The in- 
tegration measure is defined as T>a = daido2da3S(l — 
ai — 02 — 03) and the dots denote contributions of other 
Lorentz structures that drop out and also terms of twist 
5 and higher. Our notation follows Ref. [38|. 

Strictly speaking, there exist also twist-4 contributions 
from the wave function components containing two glu- 
ons or an extra quark-antiquark pair, Fig. and Fig.|3Jl, 



but they are usually assumed to be small and neglected. 
The relevant argument is based on the specific property 
of four-particle twist-4 distributions: They do not allow 
for a factorization in terms of two-particle distributions 
and, say, quark or gluon condensate. 

Higher-twist DAs can be studied using the conformal 
partial wave expansion, which is a generalization of the 
Gegenbauer polynomial expansion for the leading twist 
DAs [6l|. The contribution of the lowest conformal spin 
(asymptotic DAs) is 



50 

¥ 



Y<^(/V(i 



+ 10 )8l^)x 2 {l-xf 



xf 



(41) 



where the first and the second contribution in the paren- 
thesis are the contributions of the two-particle and three- 
particle DAs in (|38|). respectively. We stress that these 
two contributions are related by exact equations of mo- 
tion; taking into account e.g. the twist-4 correction to the 
handbag diagram and omitting contributions of three- 
particle DAs is inconsistent with QCD. One can show 
[3J] that the next-to-leading order conformal spin contri- 
butions to the relevant DAs do not contribute to ¥^(x, /i) 
so that this result is valid to NLO in the conformal expan- 
sion. The coupling 8 2 is defined by the matrix element 



(0\qgG liV 'y v q\n(p)) = ip^U&l , 

/ f,\\ 32/(9/3 ) 

«(/») = ^#rl CM. 



a s {iio) 



(42) 



9 




FIG. 3: Twist-4 corrections to the pion transition form factor 



This parameter was estimated using the QCD sum rule 
approach in [63| (see also [2fjj): 



Sl(n 2 = 1 GeV 2 ) ~ 0.2 GeV 2 . 



(43) 



D. 



Twist-six corrections 



The calculation of twist-six corrections to the transi- 
tion form factor presents a new result of this work. To 
explain why such corrections may be important, consider 
an example of the Feynman diagram shown in Fig. 2£t. 
The broken quark line with crosses stands for the quark 
condensate. If both photon virtualities are large, this 
diagram contributes to the OPE of the product of the 
two electromagnetic currents (the so-called cat-ears con- 
tribution) which involves the twist-6 four-quark pion DA 
in the factorization approximation: One quark-antiquark 
pair is put in the condensate and the other one forms a 
twist-3 chiral-odd quark-antiquark DA. Explicit calcula- 
tion gives (cf. (l7j ) 



F Fig. 4a (n 2 2\ 



where Mkillki 



V2/ w 32ira s (qq) 2 f 1 dx4>lJx) 



VfWQ 2 



xQ 2 + xq 2 



(44) 



dx ( 



(0\q(0)i>y 5 q(an)\ir(p)) = 

III,, -1- 111,1 .In 

(45) 

The DA (f^.^(x) is related to the contribution of three- 
particle (quark-antiquark-gluon) Fock state by equations 
of motion [6l| . If the contributions of the twist-3 three- 
particle DA are neglected, (f%. v {x) = 1 must be taken [H, 

EH!!. In Eq ' ® we used t ne Gell-Mann-Oakes-Renner 
relation (m u + md)((uu) + (dd)) — —f^m 2 T which is exact 
in the chiral limit. 

For the case of two equal large virtualities q 2 — Q 2 this 
contribution is of order (qq) 2 /Q 6 , so it is suppressed by 
two extra powers of 1/Q 2 compared to the leading term, 
as expected from dimension (twist) counting. On the 
other hand, the real photon limit q 2 — > of f4"4"|) cannot 
be taken in a straightforward way since the pole at q 2 = 
is clearly unphysical. This singularity appears, obviously, 
because the quark interacting with the soft photon comes 
close to the mass shell. Thus the distance it travels be- 
comes large and the OPE cannot be applied. In the full 



theory, this singularity will be tamed by nonperturba- 
tive corrections corresponding to photon emission from 
large distances, Fig. [TJd. In the simplest, vector meson 
dominance (VDM) approximation (in the LCSR method 
in addition the continuum contribution is taken into ac- 
count) nonperturbative corrections amount to a replace- 
ment 1/q 2 — > l/(m 2 + q 2 ) so that for q 2 — > the singular 
factor 1/q 2 is replaced by 1/m?. One obtains 



V2f v 2,2^a s {qqf 



dx ■ 



3 9f 2 m 2 p Q 4 J Q ~ x + xm 2 /Q 2 ' 

(46) 

Note that this is a 1/Q 4 correction to the form factor, not 
1/Q 6 as for equal virtualities. The factor 1/m 2 can be 
identified with the magnetic susceptibilit y o f the quark 
condensate (in the VDM approximation) |49l [5ol [6#j66| 



— * ~ 3.3 GcV~ 2 . 



(47) 



which enters the definition of the leading twist DA of 
a real photon [55[ so that this contribution can be 
rewritten as a convolution of photon and pion DAs with 
the coefficient function corresponding to a hard gluon 
exchange, cf. Fig. [TJd. The direct calculation of this con- 
tribution gives 



F Fi g .ib ^q2) = v^/tt 16ira s x(qq) 



3 9/2Q 4 



x / dx- 



dy- 



y 



(48) 



where 4> 7 (y) — 6y(l — y) is the leading-twist photon DA 
[49T l55j | . The integrals over the quark momentum frac- 
tions in P8")) are both logarithmically divergent at the 
end-points x — > 0, y — > 1, which signals that there is an 
overlap with the soft region, Fig. [TJ;. 

Another twist-6 contribution comes from the diagram 
in Fig. |4p 

pFi g .4b (n 2 2s V2U64Tra s (q q } 2 f 1 dx ^.Jx) 
^' W ' 9J 3 27/2 Jo (xQ 2 +xq 2 )3> 

(49) 

where the DA ^.^(x) is defined as 



•\q(0)cr^"fr o q(an)\n(p)) = ~{p^x v - p u x p ) 



m u + m d 



x / dxe-^^ix). (50) 
Jo 
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FIG. 4: Factorizable twist-6 corrections to the i ? 7 » 7 _ >w o (Q 2 ) form factor. 



In the considered approximation (neglecting the quark- 
antiquark-gluon DAs) the asy mptotic form <p^. n (x) — 
6x(l — x) must be taken [38|, [6l|, [63]. Note that this 
contribution is also of the order of l/Q 4 and not 1/Q 6 
as suggested by the naive power counting, since the limit 
q 2 — >• leads to a quadratic divergence J dx/x 2 at small 
momentum fractions. As above, this divergence is reg- 
ulated in the LCSR approach by correcting the spectral 
density to include the /9(w)-resonance and the continuum. 
Next, the diagrams in Figs. 0b, d can be calculated us- 



ing the light-cone expansion of the quark propagator in 
a background gluon field [(JT} and picking up terms con- 
taining covariant derivatives of the gluon field strength 
D^G^v. These can be reduced to a quark-antiquark pair 
via equations of motion. We have checked that there are 
no terms with additional derivatives compared to the ex- 
pression given in [67[ contributing at the required twist 6 
level. A straightforward albeit rather lengthy calculation 
gives 



l2 „2s 



V2U 167ra s (g g ) 2 f 1 f 1 u-v 

/ dv[v — v)v I du — = — 0o. (u) 
h Jo u 



27/2 
1 



1 



[uvq 2 + (1 — uv )Q 2 ] 3 [(1 — uv)q 2 + uvQ 2 ] 3 J ' 
I 



(51) 



where it was used that to our accuracy 



x (a? ( 



X (x) I m - {X) 
^ AX) 6 dx 



(52) 



Finally, the diagrams in Figs. 0£,f vanish. This is in dif- 
ference to the similar calculation for the pion electromag- 
netic form factor in Ref. [42| where only these diagrams 
contributed. 



E. The complete sum rule 

Collecting all contributions, we present here the com- 
plete light-cone sum rule with twist-6 accuracy, which 



will be used in numerical analysis in the next section. 

The sum rule for the 7r°7*7 form factor can be written 
in terms of the full QCD spectral density p(Q 2 , s) as 



V2U 



ds 



-p(Q 2 ,s) 



4 

phard 

7*7- 



so 



— 2 / dsp{Q 



2 s)eK 



-s)/M 2 



AQ 2 ) 



psoft 

7*7" 



.AQ 2 )- (53) 



Here we define the "hard" and the "soft" contributions 
as coming from large s > so and small s < so invariant 
masses in the dispersion integral, respectively. Note that 
the hard part is model-independent whereas the soft part 
is obtained under the assumption that the contribution 
of small invariant masses can be represented by a sin- 
gle narrow resonance. The continuum threshold sq can 
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be viewed as the separation scale between hard and soft 
contributions; the dependence on s has to cancel in the 
sum. 

The QCD spectral density, in turn, can be calculated 
as a sum of contributions of different twists, £=2,4,6 

p(Q 2 , s) = (Q 2 , s) + (Q 2 , a) + p^ (Q 2 , *) + ... (54) 

The leading-twist spectral density to the NLO accuracy 
is given by 



oo 



7^2 7, a n {p){ip n {x) + 



C F a s (p) 



)2 

n=0.2 



R n (Q 2 ,s)Vn(x)+i Q ^ 

fc=0,l.. 



4tt 

„ ^ H*<p k (x) 



+2 



E 

fc=0,2. 



G*<p k (x)-yy»x 



(55) 



where 

Rn(Q 2 ,s) 



= -3 



1 + 2(V>(2)- V(2 + n 

2 



(56) 



The integrals corresponding to the "hard" part, s > so, 
in Eq. ([53)) can be taken analytically. The corresponding 
expression in Eq. (E.17) in [20j contains a misprint: the 
term — 31n 2 (so/u) has to be replaced by — 31n(so/u). 
The twist-4 spectral density is equal to 

160 



P (4) (Q 2 ,*) 



3Q 



-5i(p)x i x(l-2x). 



(57) 



Finally, the twist-6 contribution can be written as: 



P (6) (Q 2 ,s) = 8i:C F a s {p) 



(qq) 2 



N c f 2 Q 6 

x + 2S(x) - - + S(x) 



2x log x + 2x log x 
1 daf 



(58) 



Note that the last two terms combine to a "plus" dis- 
tribution, 1/[1 — x] + . In all expressions (f5"5"]) -(f58 |) x = 

Q 2 /{Q 2 + s). 

We want to emphasize that the twist-6 contribution is 
not suppressed compared to the twist-4 one by an extra 
power of Q 2 , and the same is true for all higher- twist cor- 
rections. The twist expansion in LCSRs goes in powers of 
Aqcd/ s o, Aq CD /M 2 where Aq CD is a generic dimension- 
ful parameter that characterizes the size of higher-twist 
matrix elements. 



IV. NUMERICAL ANALYSIS 



326 MeV and nt = 4 active flavors. Unless stated other- 
wise, all nonperturbative parameters and models of the 
pion DA refer to the renormalization scale po = 1 GeV; 
a s (p ) = 0.494. 

A natural factorization and renormalization scale p in 
the calculation of the 7r°7*7* form factor with two large 
virtualities is given by the virtuality of the quark prop- 
agator p 2 r~j xQ 2 + xq 2 and depends on the quark mo- 
mentum fraction. If q 2 —> 0, in the LCSR framework the 
relevant factorization scale becomes p 2 ~ xQ 2 + xM 2 or 
p? ~ xQ 2 + xsq for large values of the Borel parameter, 
see e.g. [42[ ■ Note that in the first integral in (|55|) the 
quark virtuality is never large, of order Q 2 : The restric- 
tion s < so translates to x < sq/(sq + Q 2 ) and hence 
p 2 ~ 2sq as Q 2 — > oo, in agreement with the interpreta- 
tion of this term as the "soft" contribution. Numerical 
calculations with the x-dependent factorization scale are 
rather slow so in this work we use a fixed scale, replacing 
x by the constant (x) which is varied within a certain 
range: 



p 2 = (x) Q 2 + (x) s , 1/4 < (x) < 3/4 . 



(59) 



The choice of the Borel parameter in LCSRs is dis- 
cussed in [Hi], HI| • The subtlety is that the twist expan- 
sion in LCSRs goes in powers of l/(a;M 2 ) rather than 
1/M 2 in the classical SVZ approach. Hence one has to 
use somewhat larger values of M 2 compared to the QCD 
sum rules for two-point correlation functions in order to 
ensure the same hierarchy of contributions. We choose 
as the "working window" 



1 < M 2 < 2 GeV 2 



(60) 



and M 2 = 1.5 GeV 2 as the default value in our calcula- 
tions. 

We use the standard value so = 1.5 GeV 2 for the con- 
tinuum threshold as the central value, and the range 



1.3 < s < 1.7 GeV 



(61) 



in the error estimates. We did not attempt to consider 
corrections due to the finite width of the p, u> resonances. 
The estimates in Ref. [l4[ suggest that such corrections 
may result in the enhancement of the form factor by 2- 
4% in the small-to-mediate Q 2 region where the reso- 
nance part dominates. We believe that such uncertain- 
ties are effectively covered by our (conservative) choice 
of the continuum threshold. 

Finally, we use the values S 2 = 0.2 ± 0.4 GeV 2 and 
(qq) = -(240 ± 10 MeV) 3 (at the scale 1 GeV) for the 
normalization parameter for twist-4 DAs (|43[) and the 
quark condensate, respectively. 



A. The parameters 



B. Testing simple models 



All numerical results in this work are obtained us- 
ing the two-loop running QCD coupling with Aqq D = 



We start our analysis with the comparison of the LCSR 
predictions for the ir°j*j form factor for three simple 
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models of the pion DA that are often quoted in the liter- 
ature: 



6x(l — x) . 
8 



= -y/x{l - x) , 



7T 
1 



(62) 



The asymptotic 4>^(x) and flat </>^ at (a;) DAs have already 
been discussed above; the "holographic" model (j)^ ol (x) 
is inspired by the AdS/QCD correspondence |69| (see, 
however, [70l|). 

Apart from the general interest, considering these 
models allows one to test the applicability of the Gegen- 
bauer expansion. To this end, consider the approxima- 
tions to <p^ at (x), 0jj ol (x) by the truncated series at order 
n: 



^flat(hol),(n) (x) = 



flat(hol) 



fc=0,2,. 



ip k {x), (63) 



where 



hoi 



2n + 3 /T[(n + l)/2] 
3tt ^r[(n + 4)/2] 



(64) 



and a^ at are given in Eq. (fT5|) . 

The expressions collected in Sect. 3 and App. A,B al- 
low us to construct the sum rules using up to seven terms 
k = 0, 2, . . . , 12 corresponding to the n = 12 truncation. 
The resulting DAs <^ at ' (n=12) (:r), ^° l{n=12 \x) are com- 
pared with the exact ones, </>f at (x), (j)^° l (x), in Fig. [5] 
Note that the n = 12 approximation is very good for the 
"holographic" DA, whereas for the "flat" one the conver- 
gence is slow and there are large oscillations. 

Next, we use the three DAs in Eq. (1521 (the last two 
ones truncated at order n — 12) to calculate the 7r°7*7 
form factor. The results are shown in Fig. [5] One sees 
that both the asymptotic and holographic models fail to 
describe the BaBar data [l|. The flat DA fares better 
for the largest Q 2 values, but is considerably above the 
experiment at intermediate Q 2 ~ 5—15 GeV 2 . In or- 
der to understand this behavior, we compare in Fig. [7] 
the predictions for three different truncations of the fiat 
DA: 4 at > {n=12 \x), 4^ {n=8) (x) and <^ at ' (n=4) (x). We 
observe that all three calculations are very close to each 
other for Q 2 < 18 GeV 2 so that in this region contribu- 
tions of Gegenbauer polynomials starting with n = 6 — 8 
play no role. This conclusion is in agreement with our 
discussion of the MR model in Sect. 2 and also sup- 
ports the usual procedure of modelling the pion DA by 
the asymptotic expression and contributions of first two 
Gegenbauer polynomials in applications to B-decays and 
pion electromagnetic form factor, e.g. [13, fl7l - l22j . We 
also see that contributions of the Gegenbauer polyno- 
mials n = 6 , 8 become significant for the momentum 
transfers Q 2 > 18 GeV 2 . In the region Q 2 > 30 GeV 2 



0.0 0.2 



0.4 0.6 0.8 1.0 




FIG. 5: The n = 12 truncations of the pion DAs (fc° l {x) 
and 0^ at (a;) compared to the exact expressions. 




FIG. 6: The pion transition form factor for the "fiat" (solid 
red line), "holographic" (dashed blue line) and "asymptotic" 
(dash-dotted blue line) models for the pion DA, cf. Eq. (1621) . 
The experimental data are from [l|] (full circles) and 0] (open 
triangles). 



higher-order polynomials should be included into analy- 
sis as well, but the accuracy of the existing experimental 
data is not sufficient to draw definite conclusions. 

In other words, the differences between the predictions 
of asymptotic, holographic and flat DAs in Fig. [6] for 
Q 2 < 20 GeV 2 are mostly due to the different values of 
the coefficients ai and 04, with clq also playing some role. 
This leaves us with 2 — 3 parameters that can be tuned 
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FIG. 7: The pion transition form factor for three different ap- 
proximations for the "fiat" DA: n = 12 (solid red line), n — 8 
(dashed blue line),and n = 4 (dash-dotted). The experimen- 
tal data are from [lj (full circles) and (open triangles). 



Model 


scale 


«2 


04 


a 6 


a s 


aio 


ai2 


I 


fi = 1 GeV 
H = 2 GeV 


0.130 
0.089 


0.244 
0.148 


0.179 
0.097 


0.141 
0.070 


0.116 
0.054 


0.099 
0.044 


II 


/i = 1 GeV 
p = 2 GeV 


0.140 
0.096 


0.230 
0.140 


0.180 
0.098 


0.05 
0.024 


0.0 
-0.001 


0.0 

< 10" 3 


III 


jii = 1 GeV 
fi — 2 GeV 


0.160 
0.110 


0.220 
0.133 


0.080 
0.043 


0.0 
-0.001 


0.0 

< 10" 3 


0.0 

< 10" 3 



TABLE II: Gegenbauer coefficients of three sample models of 
pion DA that are consistent with BaBar measurements [1] of the 
transition form factor, cf. Fig. [9] 

to attempt a better description of the BaBar data, the 
task that we address now. 



C. Confronting the BaBar data 

The extraction of the pion DA with meaningful error 
estimates requires a global fit to the pion transition and 
electromagnetic form factor, weak B(D) — > -kIv decays 
and the couplings g^NN ,9bb*tv etc. using a Monte Carlo 
scan of the space of all available parameters, which goes 
beyond the tasks of this work. Fitting of the BaBar data 
is not attempted. Instead, we present results for three 
sample models that describe the ir°j*"f form factor suffi- 
ciently well and discuss their general features. 

The three models that we consider below are shown 
in Fig. [5] (at the scale 1 GeV) and the corresponding 
Gegenbauer coefficients are collected in Table [TTJ 

The first model 

^(a;) = l-(7/18-0.13) V a(a) (65) 

is nothing but the flat DA with the reduced second 
Gegenbauer coefficient, a% — > 0.130. It has a long 
"tail" of higher-order Gegenbauer polynomials. From 
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FIG. 8: Model I, (red solid curve), model II (blue dashed 
curve) and model III (black dash-dotted curve) of the pion 
DA at the scale 1 GeV. The asymptotic pion DA is shown by 
green dots for comparison. 



the previous discussion we expect that this "tail" actu- 
ally gives no contribution in the Q 2 range of interest. In 
order to check that this is indeed the case, we consider 
the second model in which the higher-order coefficients 
aio, a i2 are put to zero at the reference scale 1 GeV, and 
we keep a small ag to avoid an oscillating behavior at 
x ~ 1/2. Note that nonzero values of the coefficients 
Oi0i a i2 (and all higher) are generated at higher scales, 
but this mixing is numerically insignificant. Finally, the 
third model is chosen to explore the sensitivity of the 
predictions to particular values of a-2, a-4 and a^, and to 
see whether they are correlated. 

The calculations using these models are compared with 
the available experimental data in Fig. [9j The results are 
shown by thick solid curves: The line thickness shows the 
uncertainty and is calculated as a square root of the sum 
of squares of the error bars on the LCSR predictions due 
to variation of the parameters within the limits specified 
in Sect. IV. A. These include the factorization scale de- 
pendence, dependence on the Borel parameter M 2 and 
the continuum threshold srj, and on higher- twist param- 
eters 5 2 and (qq). 

The distinctive feature of all three models is the large 
value of the fourth Gegenbauer moment, 0,4, which is 
necessary in order to accommodate the observed rise of 
the scaled form factor in the Q 2 = 5 — 20 GeV 2 range. 
It is not possible to trade the large value of 04 for the 
increased 0.2 or et6, although of course there is some cor- 
relation. 

Our value of a<z ~ 0.13 — 0.16 at 1 GeV is at the low 
end of the existing estimates, cf. Table. U and in par- 
ticular it is lower compared to the earlier LCSR analysis 
of the transition form factor in Refs. [3, [20|, HI- The 
main reason for this difference is that the Borel parame- 
ter in [Tim is fixed at an ad hoc value M 2 = 0.7 GeV 2 , 
whereas in the present analysis we allow its variation in 
the 1 — 2 GeV 2 range. For the specific choice a 2 (^sy) = 
0.14, o.4(^sy) = —0.09, hsy — 2.4 GeV, advocated in 
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FIG. 9: The pion transition form factor for the three models of 
the pion DA specified in the text. The experimental data are 
from [1] (full circles) and [9] (open triangles). 



[3 Ha], the result for the form factor at Q 2 = 5 GeV 2 
is increased by ~ 11% if M 2 is changed from 0.7 to 
1.5 GeV 2 . Another reason is that in [MlMUl the twist- 
6 correction is not included. The size of this correction 
depends strongly on the Borel parameter. For our choice 
M 2 ~ 1.5 ±0.5 GeV 2 the twist-6 term proves to be small: 
factor three smaller that the twist-4 correction (see be- 
low), which is gratifying as it signals convergence of the 
OPE. In contrast, at M 2 = 0.7 GeV 2 the twist-6 correc- 
tion is almost of the same size as twist 4 and has opposite 
sign. Hence it must be included. In both cases (increas- 
ing the Borel parameter and/or including the twist-6 cor- 
rection) the net effect is the increase of the form factor by 
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FIG. 10: Contributions to the 7r°7*7 form factor from large 
( "hard" ) and small ( "soft" ) invariant masses in the dispersion 
representation, cf. Eq. (|53[) . for model I (solid curves) and 
model III (dash-dotted curves). The experimental data are 
from pi (full circles) and [f| (open triangles). 



5-10% in the CLEO range which has to be compensated 
by a smaller value of the second Gegenbauer moment. 

The error band indicated by thickness of the curves 
in Fig. [5] has to be taken with caution. A weak scale 
dependence of our results is largely due to strong can- 
cellations of the NLO radiative corrections between the 
contributions of the asymptotic DA and higher Gegen- 
bauer polynomials and may not be representative for the 
size of NNLO corrections which are only known in the 
CS factorization scheme, 



see 



for a detailed discus- 
sion of the related ambiguities. Also the uncertainty in 
the twist-4 contribution is not reduced to the (^param- 
eter: Using an alternative, renormalon model 7l| of the 
twist-4 pion DA generally produces somewhat larger cor- 
rections. We have checked that the difference is not very 
significant, however, and does not affect any of our con- 
clusions. Hence we do not show the corresponding re- 
sults. 

The "hard" and "soft" contributions to the 7r°7*7 form 
factor as defined in Eq. (f5"3")) are shown separately for 
model I (solid curves) and model III (dash-dotted curves) 
in Fig. [TUJ Asymptotically, for Q 2 — > oo, the soft con- 
tribution is power-suppressed compared to the hard one, 
~ so/Q 2 . This suppression sets in for very large values of 
Q 2 , however, especially if the pion DA is enhanced close 
to the end points. E.g. for our model III the soft con- 
tribution still accounts for ca. 25% of the form factor at 
Q 2 = 30 GeV 2 (for the separation scale s = 1.5 GeV 2 ). 
This means that a purely perturbative leading twist QCD 
calculation of the transition form factor for one real pho- 
ton in collinear factorization should not be expected to 
have high accuracy. A lattice calculation of the transi- 
tion 7rp7* form factor at Q 2 ~ 2 — 5 GeV 2 would help to 
estimate the contribution of the resonance region more 
reliably. 

Finally, in Fig. [TT] we show the higher- twist contribu- 
tions. The twist-4 correction is negative and the twist-6 
one is positive. It turns out that the twist-6 contribu- 
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FIG. 11: Higher twist contributions to the 7r°7*7 form factor 
for the values of the Borel parameter M 2 = 1.5 (solid curves) 
and M 2 = 0.7 GeV 2 (dashed curves). 
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FIG. 12: Electromagnetic pion form factor for the three models 
of pion DA described in the text. Identification of the curves fol- 
lows Fig. [8] The experimental data are from [74] (blue squares) 
and [75] (red circles). 



tion depends rather strongly on the Borel parameter. It 
is suppressed in the Q 2 region of interest relative to the 
twist-4 term for our choice M 2 =1.5 GeV 2 , but increases 
rapidly for smaller M 2 . For example, for M 2 — 0.7 GeV 2 
used in 0, [2(|, the twist-6 correction is ~ 0.6 of the 
twist-4 term at Q 2 = lGeV 2 , becomes equal (with op- 
posite sign) at Q 2 ~ 14 GeV 2 and overshoots twist-4 for 
larger Q 2 (because it contains a logarithmic ~ In Q 2 en- 
hancement). 



D. Other processes 

The pion DA is a universal function and, if extracted 
from one reaction, should, in general, describe all exclu- 
sive or semi-inclusive processes that involve a pion in the 
initial and/or final states. The most prominent of them 
are the pion electromagnetic form factor and the weak 
semileptonic decay rate B — » irivg. Without going in 
detail, we present here the corresponding LCSR calcula- 
tions using the pion DA models as specified above. 

The LCSRs for the pion electromagnetic form factor 
were derived in Refs. [H, HH, [72| and later explored also 
in 0, [zH ■ These sum rules are known to the same ac- 
curacy as for the transition form factor, i.e. including 
the NLO perturbative contribution, twist-4 and twist-6 
corrections. Explicit expressions can be found in (43[. 

The results are shown in Fig.[12]in comparison with the 
experimental data [74|, [75| • For this plot we have chosen 
M 2 = 1.5 GeV 2 , si = 0.8 GeV 2 and the factorization 
scale /i = (1/2)(Q 2 + sq) as representative values; the 
three curves correspond to the models of the pion DA in 
Fig. m 

The agreement is very good. Note that the oscillations 
at small Q 2 in model I are an artifact of the truncation 
of the Gegenbauer expansion. They can be removed, e.g. 
by reducing off* = 0.099 — > 0.04 which has the effect of 
smoothening the DA, especially in the central region. 

It has to be mentioned that the pion electromagnetic 
form factor is much more affected by soft contributions 
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Q 2 

FIG. 13: Electromagnetic pion form factor for model I (with 
reduced ai2 = 0.099 — > 0.04). The solid curve, upper dashed 
curve and lower dashed curve are calculated using the Borel pa- 
rameter M 2 = 1.5, M 2 = 2 and M 2 = 1 GeV 2 , respectively. 



The experimental data are from [7A 
circles). 



(blue squares) and [75] (red 



compared to the transition form factor and hence is also 
more model dependent. In particular the Borel parame- 
ter dependence is much stronger, see Fig. 1131 where the 
calculations using M 2 = 1 and M 2 = 2 GeV 2 are shown 
by dashed curves for comparison. 

The weak decay B — > nive has received a lot of at- 
tention as one of primary sources of information on the 
weak mixing angle \V u b\ in the Standard Model. The dif- 
ferential decay width dB/dq 2 , where q 2 is the invariant 
mass of the leptons, is given by the square of the B — > it 
form factor, modulo relevant CKM angles and kinematic 
factors 



dB 
dq 



G 2 F \V uh \ 
192tt 3 to : 



A 3/ V)|^(<7 2 )| 5 



4m|r 



(66) 
2 and 



In this equation X(q 2 ) = {m 2 B + m 2 — q 2 ) 
Tg is the mean life time of the B-meson. Below we use 
\Vub\ = 3.6 • 10~ 3 to fix the overall normalization. 

LCSRs enable one to calculate the form factor up to 




5 10 15 20 25 



Q 2 

FIG. 14: The differential -K~e + v e decay width. Experimental 
data are taken from |76j]. A 2-parameter BCL-fit [77j is applied 
to the sum rule calculation to extrapolate it to the whole q 2 
range. For normalization we assume \V u b\ = 3.6 • 1CP 3 . The 
identification of the curves follows Fig. [8] see also text. 

q 2 = 14 GeV 2 , and a two parameter BCL [77| fit is then 
used to extrapolate the calculation to the whole kine- 
matic region < q 2 < 26 A GeV 2 . The latest and most 
advanced LCSR calculations of this form factor [4(| [78| 
include NLO corrections in leading twist and also for a 
part of the twist-3 contributions. Twist-4 corrections are 
taken into account in the leading order. In the calcula- 
tions presented in Fig. Q3] we use the central values of the 
sum rule parameters from Ref. fjfjj . 

All our models of pion DA describe the data [76| 
reasonably well and are nearly indistinguishable in the 
q 2 < 14 GeV 2 range where the direct sum rule calcula- 
tion is applicable. The calculation using "conventional" 
pion DA with a 2 = 0.17 and a 4 = 0.06 (at 1 GeV) [1| is 
shown by green dots for comparison. We conclude that 
the B — > 7r form factor is not very sensitive to the higher 
Gegenbauer moments beyond a 2 ; a low value a 2 < 0.2 
(at 1 GeV) is preferred. 

The value of the pion DA in the middle point in all our 
models is close to the asymptotic value ^(x = 1/2) = 
1.5. This number is within the range in Eq. (1131) and 
somewhat larger than it was assumed in the LCSR calcu- 
lations of pion-hadron couplings (HO, l56T - l60j . The larger 
value is in fact welcome and can reduce the well-known 
~ 30% discrepancy of the sum rule calculation (5(| of 
Qd* d-k with the experiment. 



V. SUMMARY AND CONCLUSIONS 

The recent BaBar measurement [l[ of the pion transi- 
tion form factor provided one with the most direct evi- 
dence so far that the pion distribution amplitude deviates 
considerably from its asymptotic form. This result has 
to be considered as a success of an early QCD prediction 
[HI of a broad pion DA at a low scale, but it also created 
a lot of excitement because a significant scaling violation 
at Q 2 > 5 — 10 GeV 2 came out unexpected. 
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The main lesson to be learnt from the BaBar data is 
that attempts to describe the transition form factor with 
one real photon entirely in the framework of perturbativc 
QCD are futile; nonperturbative soft corrections must be 
taken into account. 

We have adopted the LCSR approach [49l - [5lT | which 
has the advantage that it is applicable to a broad class of 
reactions and has been thoroughly tested. In this work 
we go beyond the existing analysis [3, IT71422T ] in two 
aspects. First, we calculate a new, twist-six contribu- 
tion to LCSRs which proves to be sizeable. Second, we 
extend the existing formalism to allow for the contribu- 
tions of higher-order Gegenbauer polynomials, which al- 
lows one to consider DAs of arbitrary shape and also 
address the question of convergence of the Gegenbauer 
expansion which generated some confusion. 

We find that a significant rise of the scaled form factor 
Q 2 F^o rj (Q 2 ) in the Q 2 = 5-20 GeV 2 range observed by 
the BaBar collaboration |l| can be explained by a large 
value of the fourth Gegenbauer moment 

124 > (l 2 

in the pion DA, leading to models of the type shown in 
Fig. [8] which are not far from the asymptotic distribution 
in the central region but have enhancements close to the 
end points. Our preferred models also include sizeable ae 
coefficients; these can be put to zero at the cost of fur- 
ther increasing 04 which does not seem to be attractive. 
The higher partial waves, as, aio, etc. contribute only 
marginally in the BaBar Q 2 range, the reason being that 
contributions of the end-point regions in the pion DA are 
cut off by soft effects. 

We have checked that the models of pion DA having 
such an inverse hierarchy, 04 > a 2 give good description 
of the the pion electromagnetic and weak decay B — > tt 
form factors calculated within the same LCSR approach. 
This agreement is not trivial since the size of soft cor- 
rections is very different and also the duality assumption 
of the contribution of small invariant masses is applied 
in different channels. The small value of a 2 < 0.2 (at 1 
GeV) was actually suggested before from the fit to the 
B —> ■nlvt differential decay width (46|, whereas large 04 
does not have a noticeable effect in this case because the 
effective momentum transfer in B decays is much lower. 

The main uncertainty of the LCSR calculation is due to 
the assumption that contributions of low invariant masses 
in the dispersion relation in QCD diagrams are dual (i.e. 
coincide in integral sense) with the contribution of reso- 
nances, here the p, uj mesons. The accuracy of duality is 
difficult to quantify, but it is usually believed to be bet- 
ter than 20% on the experience of many successful ap- 
plications. An inspection shows that our result 04 > a 2 
is related to a rather large value of the transition form 
factor F pvi *(Q 2 ) in the Q 2 ~ 2 — 5 GeV 2 range that 
follows from duality. For comparison, this form factor 
estimated as the integral of the spectral density below 
s = 1.5 GeV 2 in the MR model [23[ appears to be a 
factor 2-3 lower, which explains why in this model the 
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BaBar data can be fitted by a CZ-type pion DA with 
a large a 2 coefficient and = 0. Lattice calculations 
of the Fpiry* (Q 2 ) form factor in a few GeV 2 range and 
improved accuracy on 02 would help to discriminate be- 
tween these two possibilities. More precise experimental 
data in the Q 2 = 15 — 30 GeV 2 range would of course be 
most welcome as well. 
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The NLO anomalous dimensions can most easily be 
obtained using the FeynCalc Mathematica package [79j . 
For convenience we present explicit expressions up to n = 
12 that are used in our calculations (7Q 1 ' = 0): 



7^ 
7$ 



34450 830 



243 8l" n/ ' 
662846 31132 



-n/, 



3375 2025 
718751707 3745727 



-n f , 



3087000 198450 
293323294583 19247947 



-n f , 



1125211500 893025 
212204133652373 512808781 



-n f , 



748828253250 21611205 
995653107122188087 93360116539 



3290351344780500 



3652293645 



n f . 



The off-diagonal mixing coefficients dfj in Eq. ^ are 
given by the following expression: 



Appendix A: Scale dependence of the pion DA 

The scale dependence of the coefficients a n (fj) in the 
Gegenbauer expansion of the pion DA is determined by 
Eq. ©. 

The RG factor E^ LO (^i, jUq) in this expression is given 

by 



7i 0) /2/3o 



X 1 + 



a s (Mo)_ 
a s (ii) - a s (no) 7i 0) / 7^ 



(A.l) 



8tt 



A) V7i 0) A) 




The corresponding LO RG factor is obtained 

by keeping the first term only in the braces. 

Here ft (ft) and 7^ (7^ ') are the LO (NLO) coeffi- 
cients of the QCD ftfunction and the anomalous dimen- 
sions, respectively: 



1 da a (iJ.) 
dfi 2 



fta s ) = -a s {^+ft(^) 2 + -} 



ln[as) ~ 2l 7n 4n +ln UJ 



The first two coefficients of the beta-function are 

2 38 
ft, = ll-~ri/, ft = 102-yn/, 

whereas jn is given by 



n+l 



In 



i 0) = 2C F 1 - 



in + l)(n + 2) — - in 



4]T- 

' J 711 



(A.2) 



(A.3) 



(A.4) 



d*(/i,Mo) 



7n 



1 



(0) 

a s (Mo) 



2ft 

bi° ) -7f ) -2^o]/2/3 



(A.5) 



The matrix Af^ is defined as 
, fc _ (fc + l)(fc + 2)(2n + 3) 



(n + l)(n + 2) 
(0) - 2/3 



1 (n- k)in + k + 3) 



Ad 



^-^(n + 2) + ^(l) 



(fc + !)(* + 2) 



(A, 



where 



A, 



+ 2i/)(n-k)-ip(n + 2)-ip(l). (A.7) 



For convenience, we have collected numerical values of 
the coefficients M% for n < 12 in Table [En] 



Appendix B: The NLO spectral density 

The coefficients G\ and if„ in the expansion of the 
NLO perturbative spectral density ([50)1 are collected in 
Tabs. [IV] and E] respectively. Our results for G^ agree 
with Ref. [3 (except for G% and G°) noting an overall 
sign difference in definition of G\ , whereas for H% the 
difference is that the expansion in Eq. ([50)1 also involves 
contributions with odd k = 21 + 1. 
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fc = 


k = 2 


k = 4 


fc = 6 


k = 8 


fc = 10 


n = 















77 = 2 


-11.23 4- 1.7371/ 













n = 4 


-I.4l4-0.56n/ 


-22.02 + 1.65n/ 











n = 6 


0.03 4-0.26?!/ 


-7.76 4- 0.8277/ 


-22.77 4- 1.3977/ 









n = 8 


0.29 4-0.14n/ 


-3.34 + 0.4877/ 


-10.34 + 0.8477/ 


-21.72 + 1.1877/ 







n= 10 


0.31 4-0.0971/ 


-1.58 4- 0.3077/ 


-5.46 + 0.5577/ 


-11.3 + 0.7977/ 


-20.35 + 1.0277/ 





n = 12 


0.28 4-0.0677/ 


-0.78 4- 0.2177/ 


-3.13 + 0.3877/ 


-6.64 + 0.5677/ 


-11.54 + 0.7377/ 


-19.0 + 0.977/ 



TABLE III: The mixing matrix M„ flAp . 





fc = 


k = 2 


k = 4 


k = 6 


fc = 8 


fc = 10 


fc = 12 


77 = 


-1 














77 = 2 


3 
2 


35 
12 












77 = 4 


3 
4 


161 

72 


203 
45 










77 = 6 


83 
180 


49 
40 


781 
300 


29531 
5040 








n = 8 


177 
560 


4 
5 


6259 
4200 


4437 
1568 


177133 
25200 






77 = 10 


487 
2100 


6181 
10800 


7601 
7560 


7823 
4704 


338561 
113400 


1676701 
207900 




77 = 12 


74141 
415800 


17167 
39600 


697 
945 


177799 
155232 


2227921 
1247400 


5672237 
1829520 


30946717 
3439800 



TABLE IV: Numerical values of the coefficients G„ <[33|). 
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